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Abstract Comparison of appropriate theoretical derived line ratios with observational 
data can yield estimates of a plasma's physical parameters, such as electron density or 
temperature. The usual practice in the calculation of the line ratio is the assumption of 
excitation by electrons/protons followed by radiative decay. Furthermore, it is normal to 
use the so-called coronal approximation, i.e. one only considers ionization and recom- 
bination to and from the ground-state. A more accurate treatment is to include the ion- 
ization/recombination to and from meta-stable levels. Here, we apply this to two lines 
from adjacent ionization stages; Mg IX 368 A and Mg X 625 A, which has been shown 
to be a very useful temperature diagnostic. At densities typical of coronal hole condi- 
tions, the difference between the electron temperature derived assuming the zero density 
limit compared with the electron density dependent ionization/recombination is small. 
This however is not the case for flares where the electron density is orders of magnitude 
larger. The derived temperature for the coronal hole at solar maximum is around 1 .04 MK 
compared to just below 0.82 MK at solar minimum. 
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1 INTRODUCTION 

Over the past few decades, many electron temperature and density diagnostics have been developed for 
application to observations of lines in the ultraviolet and extreme-ultraviolet regions of the spectrum, 
thereby providing insight to the physical parameters of the solar transition region and corona. The line 
ratio method usually involves a pair of lines from the same ion, although lines from adjacent ionization 
stages may be used. Thus, provided we have an accurate atomic model, the derived temperature provides 
essential data for an understanding of the basic processes occurring on the Sun. 

Coronal holes, seen as regions with reduced emission due to a lower electron temperature than 
the surrounding quiet Sun, are an important source of open magnetic flux and high-speed solar wind 
plasma. Furthermore, observing the evolution of coronal holes over the course of a solar cycle can lead 
to a better understanding of the solar dynamo and advance our knowledge in forecasting space weather. 
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However, all this requires an effective method to detect and track coronal hole boundaries. Harvey & 
Recely (2002) used a hand-drawn method involving the He I 10830 A line plus photospheric magnetic 
field data. This was further developed by Henney & Harvey (2005) and Malanushenko & Jones (2005) 
showing good agreement between the hand drawn method and SoHO/EIT 195 A images. With the 
amount and quality of EUV data currently available from SoHO, a few attempts have been made to 
develop a more automated method of coronal hole identification. 

Scholl & Habbal (2008) used SoHO/EIT Fe ix/x 171 A, Fe XII 195 A and He II 304 A, 
plus SoHO/MDI magnetograms. Chapman (2008, 2009) used EUV spectral line data obtained with 
SoHO/CDS. He showed that the ratio of Mg IX 368 A and Mg X 625 A resonance lines are an excellent 
temperature diagnostic for coronal holes, thus enabling a method not only to locate the coronal hole 
boundary but also an evaluation of the electron temperature with time. This work however employed 
the so-called coronal approximation, or zero-density approximation. This treats the populations of ex- 
cited states of ions via an excitation balance of collisional excitation by electrons, and radiative decay. 
The ionization state is established as a balance of electron impact ionization from the ground state and 
radiative plus dielectronic recombination. In the simplest version of such modelling, secondary colli- 
sions with excited states are neglected. Following on from work by Burgess & Summers (1969) and 
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Fig. 1 The observed ratio of Mg IX 368 A to Mg X 625 A from 1998 onwards as derived by 
Chapman (2008, 2009) smoothed over 28 days. The typical error on this smoothed data is less 
than 2%. 



Vernazza & Raymond (1979), it was shown by Doyle et al. (2005, 2006) that the contribution functions 
of resonance lines from low Z Be-like and Li-like ions can be altered substantially due to electron den- 
sity dependent ionization considerations. Although, we know that coronal holes are regions of reduced 
electron density compared to the quiet Sun, we do not know the extent of this additional atomic process 
with regard to Mg IX and Mg X. 

Here, we investigate this dependence for Mg and apply it to the synoptic data obtained with 
SoHO/CDS (Harrison et al. 1995) as discussed by Chapman (2008, 2009) in order to access the ex- 
tent of the electron density sensitivity in the Mg IX 368 A and Mg X 625 A lines. 



A coronal temperature diagnostic 



3 



2 CONTRIBUTION FUNCTION 

An important component in calculating the line contribution function is the ionization balance. As would 
be expected of analyses of the solar corona, the so-called coronal approximation is most often imple- 
mented when calculating the ionization balance of the emitting plasma. In this approximation the plasma 
is considered optically-thin, low-density, dust-free, and in steady-state or quasi-steady-state. Under these 
conditions the effects of any radiation field can be ignored, three-body collisions are unimportant, and 
the ionization balance of the gas is time-independent. The population of ionization stages is calculated 
as a balance of ionization due to electron impact ionization (EII) and recombination due to radiative re- 
combination (RR) and dielectronic recombination (DR). The most recent compilation of the ionization 
balance under these conditions was by Bryans et al. (2009). 

However, under certain conditions, the low-density approximation may not be appropriate as activ- 
ity in the solar atmosphere can result in an increase in density where the coronal approximation breaks 
down. In the coronal approximation there is no treatment of metastable states with populations compa- 
rable to the ground state. In addition, the coronal approximation takes no account of secondary electron 
collisions. The generalized collisional-radiative picture (Summers & Hooper, 1983) allows such an anal- 
ysis. In detail, the collisional ionization and redistribution processes from excited states are included. 
The populated metastable states are determined via an elaborated ionization balance along with the 
ground states. These were computed within the Atomic Data and Analysis Structure (ADAS; Summers 
20093) framework, which is a collection of fundamental and derived atomic data, and codes that manip- 
ulate them. 

Considering metastable level populations (including the ground level) to be dynamic and excited 
levels (i.e. all other levels) to be quasi-static, one can form population equations for each N p and 
Ni'in terms of these matrix elements as 

C pa N^ + C PJ Nf + C pa+ N^ 1) + C pa -N { ;- 1] (1) 

C ia N^ + dj + C ia+ Nlf 1} + C iCT _ (2) 

where we denote metastable states (including the ground state) by Greek indices and excited states by 
Roman indices. In the above we use the summation convention over repeated indices. Such a matrix 
contains elements C pa , C p j, Ci a , and Cy denoting the excitation (and de-excitation) rates between the 
respective levels; C pa+ and Ci a+ , denoting the sum of the recombination rates from each metastable; 
and C pa _ and Gi a _, denoting ionization rates; while the diagonal elements C pp and Cu, indicate total 
collisional and radiative loss rates from the levels p and i. 

Recent work by Doyle et al. (2005, 2006) looked at this problem regarding formation of the res- 
onance lines from low Z Li-like and Be-like ions. In the present work we use ADAS to calculate the 
contribution functions of Mg IX 368 A and Mg X 625 A for different values of the electron density. The 
results of the ratio of these contribution functions are shown in Fig. 2. We use the ionization and recom- 
bination data from ADAS in these calculations. The variation of the contribution functions is primarily 
due to the suppression of dielectronic recombination with increasing density as collisional depopulation 
of high n-shells, formed after stabilization, takes effect. This results in the ionization balance being 
shifted to lower temperatures and is reflected in the contribution function. In the next section, we apply 
the above results to recent coronal hole data. 

3 OBSERVATIONAL DATA 

The observational data used by Chapman and reproduced here (see Fig. 1) were the CDS synoptic study 
(SYNOP) which created 2 dimensional images and were run daily from around mid-night. The study 
generated a series of nine rasters each of which covers a 4x4 arcmin area. Together, the nine rasters 
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Fig. 2 The ratio of the electron density dependent contribution function for Mg IX 368 A 
and Mg X 625 A for four values of electron density ranging from very low densities to flare- 
type densities. The short horizontal lines at ratios of 1.0, 2.5 and 4.0 are the observed solar 
minimum, mean and maximum values, respectively (see Fig. 1). 



spanned the central meridian of the Sun. The study used a selection of the strongest lines available in 
the CDS/N1S wavelength range, covering a wide range of temperatures from the chromosphere to the 
corona: He I 584 A, O V 629 A, Mg VII 367 A, Mg IX 368 A, Mg X 625 A and Fe XVI 361 A. As a 
result, a near-continuous dataset exists since 1996 which is well-suited to studying long-term variations 
in the solar structure. However, here we only use data taken after the SoHO recovery, i.e. from 1998; 
due to cross-calibration problems which are discussed in more detail by Chapman (2008, 2009). Various 
methods used to extract the line fluxes, including background subtraction, use of small windows, line 
fitting, calibration, are all discussed by Chapman (2008, 2009). 

4 RESULTS 

At densities typical of coronal hole conditions (i.e. 10 8 cm -3 to 10 9 cm~ 3 , see Raymond & Doyle 1981), 
the difference of including ionization/recombination to/from meta-stable levels is small compared to the 
normal coronal approximation. This however, is not the case for flares where the electron density can 
be orders of magnitude larger. In Fig. 3, we show the derived electron temperature over a 9 year period 
from 1998 onwards for different electron densities. Note the lower temperature at solar minimum (w 
2008), while at solar maximum (« 2001) there are intervals where the electron temperature is above 
1 MK. Due to the high quality of the observational data, these observed changes, including not only 
the difference at solar minimum/maximum, but also the short-term changes throughout the cycle are 
real. The paper by Chapman (2009) discusses in more detail the changes in the coronal hole coverage 
throughout the solar cycle, and its implications. 
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Fig. 3 The derived electron temperature over a 9 year period from 1998 onwards for different 
electron densities. Note the lower temperature at solar minimum (w 2007/8), while at solar 
maximum (w 2001/2) there are intervals where the electron temperature is above 1 MK. 
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